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In this work we address the recent experiments of Altimiras and collaborators,—^ where an electron 
distribution function at the quantum Hall (QH) edge at filling factor v — 2 has been measured with 
high precision. It has been reported that the energy of electrons injected into one of the two chiral 
edge channels with the help of a quantum point contact (QPC) is equally distributed between them, 
in agreement with earlier predictions, one being based on the Fermi gas approach^ and the other 
utilizing the Luttinger liquid theory.— We argue that the physics of the energy relaxation process 
at the QH edge may in fact be more rich, providing the possibility for discriminating between two 
physical pictures in experiment. Namely, using the recently proposed non-equilibrium bosonization 
technique^ we evaluate the electron distribution function and find that the initial "double-step" 
distribution created at a QPC evolves through several intermediate asymptotics, before reaching 
eventual equilibrium state. At short distances the distribution function is found to be asymmetric 
due to non-Gaussian current noise effects. At larger distances, where noise becomes Gaussian, the 
distribution function acquires symmetric Lorentzian shape. Importantly, in the regime of low QPC 
transparencies T the width of the Lorentzian scales linearly with T, in contrast to the case of equi- 
librium Fermi distribution, whose width scales as y/T. Therefore, we propose to do measurements 
at low QPC transparencies. We suggest that the missing energy parados^ may be explained by the 
non-linearities in the spectrum of edge states. 

PACS numbers: 73.23.-b, 03.65.Yz, 85.35.Ds 



I. INTRODUCTION 

A two dimensional electron gas (2DEG) in strong per- 
pendicular magnetic field exhibits the regime of quan- 
tum Hall effect^ (QHE). One of the peculiar phenom- 
ena specific to this regime is the appearance of one di- 
mensional (ID) chiral edge states, which are quantum 
analogs of skipping orbits. Recent extensive experimen- 
tal studies^— of these sates have led to the emergence of 
a new field in condensed matter physics dubbed the elec- 
tron optics. On the theoretical side, there are two main 
points of view on the physics of quantum Hall (QH) edge 
states. One group of theories^ suggests that at integer 
values of the Landau levels filling factor the edge exci- 
tations are free chiral fermions. The second group of 
theories is based on the concept of the edge magneto- 
plasmon picture.™ The fundamental edge excitations in 
these theories are the charged and neutral collective bo- 
son modes. 

The domain where these two approaches meet each 
other is the low-energy effective theory?^ In the frame- 
work of this theory, both fermion and boson excitations 
are two forms of the same entity. Namely, they can be 
equivalently rewritten in terms of each other: 

ip(x,t) ~ exp[i(p(x,t)] 

where ijj(x,t) is the fermion field, and </>(x, i) is the bo- 
son field. However, this transformation is highly nonlin- 
ear, and in the presence of strong Coulomb interaction 
fermions are not stable and decay into the boson modes 
which arc the cigcnstatcs of the edge Hamiltonian. 

Results of tunneling spectroscopy experiments^ rea- 
sonably agree with the free-electron description of edge 



states. However, the first experiment on Aharonov-Bohm 
(AB) oscillations of a current through the electronic 
Mach-Zehnder (MZ) interferometer— has shown that the 
phase coherence of edge states is strongly suppressed at 
energies, which are inversely proportional to the interfer- 
ometer's size. Moreover, several subsequent experiments 
on MZ interferometers at filling factor v = 2 have shown 
puzzling results on finite bias dephasingSr— theoretically 
studied ir>i£~— . Namely, the visibility of AB oscillations 
in these experiments is found to have a lobe-type pat- 
tern as a function of the applied voltage bias. Such re- 
sults are difficult to explain in terms of the fermion pic- 
ture, while they all follow naturally from the plasmon 
physics, 1? where the Coulomb interaction plays a crucial 
role. Thus the boson picture of edge excitations might 
be more appropriate. 

In contrast to mentioned above non-local experiments, 
some local measurements seem to be not able to differ- 
entiate between two physical pictures of edge states. For 
example, both theories predict Ohmic behavior of the 
tunneling current, unless it is renormalized by a non- 
linear dispersion of plasmons. Moreover, the equilibrium 
distribution of the bosons is equivalent to the one of 
fermions (see the demonstration of this fact in Section 
IIV Cj) . Therefore, it might be interesting to investigate 
non-equilibrium local properties of edge states. 

Non-equilibrium behavior of ID systems has been 
a subject of intensive theoretical 2 ^ and experimental 22 
studies for a long time. However, only recently it has 
become possible to measure an electron distribution at 
quantum Hall edge /(e) as a function of energy e with 
high precision.™ The main idea of the experimental tech- 
nique is to restore the function /(e) by measuring the dif- 
ferential conductance Q of tunneling between two edges 
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through a single level in a quantum dot: 
G(e) cx df(e)/de, 



(1) 



where e is the energy of the quantum dot level, controlled 
by the gate voltage V g . This technique has been used in 
experiments Q in order to investigate the energy relax- 
ation at QH edge states at filling factor v = 2. The 
schematics of these experiments is shown in Fig. [TJ The 
main result is that the electron distribution relaxes to- 
ward local equilibrium Fermi distribution, and the energy 
splits equally between the two edge channels. 




FIG. 1: Schematics of the experiments pj and Q. The shaded 
region is filled by the 2D electron gas in the regime of the 
quantum Hall effect. At filling factor v = 2 there are two chi- 
ral edge states shown by the blue (the outer channel) and the 
black (the inner channel) lines. The QPC of the transparency 
T and biased with the voltage difference A/x, injects electrons 
into the outer channel, and thus creates a non-equilibrium 
electron distribution. After the propagation along the QH 
edge, the distribution is detected at distances L from the 
source with the help of a quantum dot with a single level 
controlled by the gate voltage V g . 

The first theoretical models, based on the fermion pic- 
ture Q and on the plasmon approach Q , have come qual- 
itatively to identical conclusions. Namely, both works 
predict equal distribution of the energy between the edge 
channels, in agreement with the experimental findings. 
In other words, based on the results of Refs. Q and Q 
alone, the experimentalists are not able^ to discriminate 
between two alternative descriptions of the physics of QH 
edge. Thus, it seems to be important and timely to re- 
analyze the problem of the energy relaxation at the QH 
edge in order to make new, model specific and distinct 
predictions that can be verified experimentally. This is 
exactly the purpose of the present work. 

Here we show that the Coulomb interaction strongly 
affects the spectrum of collective edge excitations and 
leads to the formation of charged and dipole plasmons 
modes, which propagate with different velocities^ They 
carry away the energy of electrons injected through the 
QPC and equally distribute it between edge channels at 
distances L ex from the QPC. In addition to this observa- 
tion, which agrees with findings of previous works^i! we 



stress that the same process splits the wave packets of 
injected electrons, and leads to strong coupling of elec- 
trons to the noise of the QPC current. The regime of 
weak injection, i.e., when the transparency of the QPC 
is small, TCl, deserves a special consideration. In this 
regime the current noise at relevant time scales becomes 
markovian, and as a result, the function —df(e)/de ac- 
quires a Lorentzian shape. (This effect resembles a well 
known phenomenon of the homogeneous level broaden- 
ing.) Interestingly, the width of the Lorentzian scales as 
TA/i at small T, where A/i is the voltage bias applied 
to a QPC. In contrast, the width of the eventual equilib- 
rium Fermi distribution of thermalized electrons scales as 
s/TA/j,. If thermalization takes place at longer distances, 
L cq 3> L C x, then the intermediate regime described here 
may be observed in experiment with a weak injection. 
This would indicate that interactions strongly affects the 
physics at the edge and that the fermion picture becomes 
inappropriate. 

In order to theoretically describe the experiments Q 
and to quantitatively elaborate the physical picture, we 
use the non- equilibrium bosonization technique, which 
has been introduced in our previous work The main 
idea of this approach is based on the fact that in a ID 
chiral system one can find a non-equilibrium density ma- 
trix by solving equations of motion for plasmons with 
non-trivial boundary conditions. Then, one can rewrite 
an average over the non-equilibrium state of an interact- 
ing system in terms of the full counting statistics (FCS) 
generators^! of the current at the boundary. In the sit- 
uation considered in this paper, because of chirality of 
QH edge states, interactions do not affect the transport 
through the QPC alone. This leads to a great simplifica- 
tion, because in the markovian limit the FCS generator 
for free electrons is known^l 

The structure of the paper is following: In Sec. [XT] we 
describe the non-equilibrium bosonization technique in 
some details. Next, we use this technique in Sec. Hill in 
order to find the electron correlation function for different 
distances from the QPC. Finally, we use these results to 
find the electron distribution function in Sec. IIV1 and 
present our conclusions in Sec. IVII Several important 
technical steps and the phenomena resulting from the 
non-linearity of the spectrum of plasmons are described 
in Sec. IVland Appendices. 



II. NON-EQUILIBRIUM BOSONIZATION 



We note, that the relevant energy scales in recent meso- 
scopic experiments with the QH edge statei£~— are very 
small. Therefore, it is appropriate to use the low-energy 
effective theory^ of the QH edge. One of the advantages 
of this theory is that it allows to take into account strong 
Coulomb interactions in a straightforward way^ How- 
ever, an additional complication arises from the fact that 
in experiments^ the injection into one of the two edge 
channels creates a strongly non-equilibrium state. We, 
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therefore, start by recalling in this section the method 
of non-equilibrium bosonization, proposed earlier in Ref. 

which is suitable for solving the type of a problem 
that we face. Throughout the paper, we set e = h = 1. 



A. Fields and Hamiltonian 

According to the effective theory of QH edge^ the 
collective fluctuations of the charge densities p a (x) of the 
two edge channels, a = 1, 2, at filling factor v = 2 are the 
only relevant degrees of freedom at low energies. These 
charge densities may be expressed in terms of the chiral 
boson fields, <fi a (x), 

p a (x) = (l/2n)d m (l> a (x), (2) 

which satisfy the following commutation relations: 

[<t> a (x) , (/)/} (y)} = iir5 a psgn(x - y) , (3) 

The vertex operator 

1 



(4) 



annihilates an electron at point x in the edge channel a. 
The constant a in the prefactor is the ultraviolet cutoff, 
which is not universal and will be omitted and replaced 
by other normalizations. One can easily check, with the 
help of the commutation relations Q, that the operators 
@ indeed create a local charge of the value 1 at point x, 
and satisfy fermionic commutation relations. 

Close to the Fermi level the spectrum of electrons may 
be linearized, therefore the free-fermion part Ho of the 
total QH edge Hamiltonian H = Hq + Hint takes the 
following form: 



Hq = —ivp 



dxipl,d x ^ a , 



(5) 



where the bare Fermi velocity vf is assumed to be the 
same for electrons at both edge channels. The sec- 
ond contribution to the edge Hamiltonian describes the 
density-density Coulomb interaction, 

H m t = (1/2) ^2 dxdy U <*p( x " y)Pa(x)Pfs(y), (6) 

which is assumed to be screened at distances d smaller 
than the characteristic length scale L in experiments [7J, 
HHlll, i.e., L > d. Therefore, we may write: 



U a p{x -y) = U a fi5(x - y). 



(7) 



Screening may occur due to the presence of either a back 
gate, or several top gates. We show below that the as- 
sumption (J7J results in the linear spectrum of charge ex- 
citations. This approximation seems to be reasonable, 
agrees well with some experimental observations such as 



an Ohmic behavior of the QPC conductance at low volt- 
age bias, and eventually does not strongly affect our main 
results. Nevertheless, below we relax this assumption and 
investigate the effects of weak and strong non-linearities 
in the spectrum of charge excitations. 

After taking into account the relations © and (J3J) and 
applying the point splitting procedure, we arrive at the 
edge Hamiltonian of the quadratic form in boson fields: 



H = g^y^ Vqfs ( dxd x <p a {x)d x ct)i 3 {x) 

a,B J 



(8) 



which nevertheless contains free fermion contribution as 
well as the Coulomb interaction potential: 



V a /3 = 2TTV F S a p + U a p. 



(9) 



Equations §5§, (|S]) and ([!} complete the description 
of the QH edge at low energies. 

The experimentally founds electron distribution func- 
tion at the outmost QH edge channel is given by the ex- 
pression: 

/(e) = [dte~ iet <v4(£,t)lML,0)). (10) 



Rewriting this expression via the boson fields we finally 
obtain 



f(e) = J dte- M K{t) 
K{t) = 



i0i(L,t) e i(^i(L,O)\ 



(11a) 
(lib) 



where we have introduced the electron correlation func- 
tion K, evaluated at coincident points at distance L from 
the QPC. The proportionality coefficient in (lllb[) may be 
corrected later from the condition that /(e) takes a value 
1 for energies well below the Fermi level (see, however, the 
discussion in Sec. [V] for further details). In equilibrium, 
in order to evaluate the correlation function on the right 
hand side of this equation one may now follow a standard 
procedure^ of imposing periodic boundary conditions on 
the boson fields and diagonalizing the Hamiltonian ©. 
In our case, however, the average in (|11[) has to be taken 
over a non- equilibrium state created by a QPC. Attempt- 
ing to express such a state entirely in terms of bosonic 
degrees of freedom is a complicated and not a best way 
to proceed. We circumvent this difficulty by applying a 
non-equilibrium bosonization technique proposed in our 
earlier work _5] . This technique is outlined below in some 
detail. 



B. Equations of motion, boundary conditions, 
and FCS 

The Hamiltonian ([8]), together with the commuta- 
tion relations ([3]), generates equations of motion for the 
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fields <fi a , which have to be complemented with boundary 
conditions^ 

d t (j> a (x,t) = -—^2V a0 d x <pp(x,t), (12a) 

d t <f> a (0,t) = -2nj a (t). (12b) 

The last equation follows from the charge continuity con- 
dition d t p a + d x j a — and the definition ([2]). Thus the 
operator j a (t) describes a current through the bound- 
ary x = in the channel a. For the convenience, we 
place a QPC in the outer channel a = 1 right before the 
boundary, so that the operator j\ (t) describes an outgo- 
ing QPC's current. 

The key idea of the non-equilibrium bosonization ap- 
proach is to replace the average in Eq. (ITT]) by the average 
over temporal fluctuations of currents j a , the statistics of 
which is assumed to be known. Indeed, although in gen- 
eral the fields <f> a influence fluctuations of the currents 
j a , leading to such effects as the dynamical Coulomb 
blockade^! and cascade corrections to noise,— in the case 
of chiral fields describing QH edge states no back-action 
effects arise i 16 ' 18 As a consequence, at integer filling fac- 
tors the electron transport through a single QPC is not 
affected by interactions, which seem to be an experimen- 
tal factjiiii Therefore, by solving equations (IT21 one may 
express the correlation functions of the fields 4> a in terms 
of the generator of full counting statistics (FCS): 24 



_ / e »AQ Q (t) e - l AQ ct (0)\ 



(13) 



Here averaging is taken over free electrons, and the op- 
erators 



Qa(t) 



dt'j a (t' 



(14) 



may be viewed as a total charge in the channel a to the 
right of the boundary at x = 0. 

To prove the connection of the electron correlations 
in (fTTj) to the generating functions (|13p . we come back 
to the discussion of the interaction effects, which are in 
fact encoded in a solution of the equations of motion 
(I12ap . The long-range character of the Coulomb interac- 
tion leads to the logarithmic dispersion in the spectrum 
of collective charge excitations, the physical consequences 
of which are discussed in Sees. HVI andlVl For a moment, 
to simplify equations (I12a|l . we have assumed screening 
of the Coulomb potential at distances d shorter than the 
characteristic length scale Vp/Afj,, which is of the order 
of few microns in recent experiments. Nevertheless, it 
is very natural to assume that the screening length d is 
much larger than the distance a between edge channels, 
d^$> a, which does not exceed few hundreds nanometers. 
Therefore, one can write 



U a p = vru, u/vp ~ log(d/a) > 1, 



(15) 



i.e., the in-channel interaction strength is approximately 
equal to the intra-channel. As a result, the spectrum of 



collective charge excitations splits into two modes: a fast 
charged mode <fii with the speed u, and a slow dipole 
mode 4>2 with the speed v ~ vp. 

It is important to stress that the condition c!> a, lead- 
ing to (|15p , results in a sort of universality: the solution 
of equations of motion (|12ap in terms of the charged and 
dipole mode, 

4>i (x, t) = —= [0i (x — ut) + <f> 2 (x — vt)] (16a) 
v2 

4> 2 (x, t) = — [(j>i(x - ut) - <f> 2 (x - vt)] (16b) 
v2 

is only weakly sensitive to perturbations of our model, in 
particular to those that account for different bare Fermi 
velocities of edge channels and slightly different interac- 
tion strengths. 

Applying now boundary conditions (|12bp to the result 
(|16l) , we finally solve equations of motion in terms of the 
boundary currents: 



<MM) = -W dt'\j 1 (t')+j 2 {t')] 



W dt'[ 3l (t')-j 2 (t% (17a) 



<t> 2 {x,t) = -n dt'\ji(t') + h{t')] 



where we have introduced notations 

t u = t — x/u. t v — t — x/v. 



+ n dtW) - j 2 (t% (17b) 



(18) 



Finally, using the definition (1141) , we arrive at the solution 
in the compact form 

<t>i(x,t) = -n[Ql(tu) + Q2(t u ) + Ql(tv)-Q2(tv)}, (19) 

and to a similar expression for the inner channel. The 
physical meaning of this result is rather simple: when 
charges are injected into the channel a — 1 and 2, they 
excite charged and dipole mode (note the minus sign in 
the fourth term on the right hand side) which have dif- 
ferent propagation speeds u and v. As a result, these 
charges arrive at the observation point x with different 
time delays x/u and x/v, and make a contribution to the 
field 4>\ at different times (fT5|) . 

When substituting this result into the correlation func- 
tion in Eq. (jllbp one may use the statistical independence 
of the current fluctuations at different channels and split 
the exponential functions accordingly: 



K(t) 



/^[Qii,t u )+Qi{t v )] e -in[Qi(t u -t)+Qi(t v -t)]\ 
x / e iTr[Q2(t u )-Q 2 (t v )] e -i-K[Q 2 (t v .-t)-Q 2 {t v -t)]\ (20) 



In the rest of the paper we will be interested in the cor- 
relation function at relatively long distances L 3> vt c , 
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FIG. 2: Schematic illustration of the Coulomb interaction ef- 
fect at the HQ edge at filling factor v — 2. The electron 
wave-packet of the charge e created in the outer edge channel 
(black line) decays into two eigenmodes of the Hamiltonian 
p|). the charged and dipole mode, which propagate with dif- 
ferent speeds and carry the charge e/2 in the outer channel. 
As a result, the wave packets do not overlap at distances larger 
than their width, and contribute independently to the elec- 
tron correlation function with the coupling constant A = it. 29 
Similar situation arises when an electron is injected in the 
inner channel (blue line), however in this case the charged 
and dipole states carry opposite charges at the outer chan- 
nel. Thus, there are four independent contributions to the 
correlation function in the outer edge channel. 

where t c ~ 1/A/i is the correlation time of fluctuations 
of the current through a QPC. (We show below that at 
this length scale the energy exchange between two chan- 
nels takes place.) In this case, the partitioned charges 
Q a , taken at different times t u and t v , are approximately 
not correlated, as illustrated in Fig. [3] This assumption 
is quite intuitive and may be easily checked using Gaus- 
sian approximation. We, finally, arrive at the following 
important result: 



tf(t)=Xi(M)X2(-7r,t)x2(7r,t), 



(21) 



i.e., the electronic correlation function (|20p may indeed 
be expressed in terms of the FCS generator (fTS)) , 



III. ELECTRON CORRELATION FUNCTION 

The expression (j2"Tj) presents formally a full solution 
of the problem of evaluation of an electron correlation 
function. Generators of the FCS for free electrons in 
this expression, defined as (|13|) . may be represented as 
a determinant of a single particle operator^ and even- 
tually, evaluated, e.g., numerically. However, a further 
analytical progress is possible in a number of situations, 
which are important for understanding physics of the en- 
ergy relaxation processes. In particular, we show in this 
section that for the case of equilibrium fluctuations of 
the boundary currents, the correlation function ([2~Tj) as 
well as the electron distribution function (|11[) acquire an 
equilibrium free-fermionic form. The electron correla- 
tion function may also be found analytically away from 
equilibrium for the case of a Gaussian noise. Interest- 
ingly, in the short-time limit, t <C 1/A/i, the main con- 
tribution to the correlation function comes from zero- 
point fluctuations of boundary currents, and it behaves 
as a free-fermion correlator, i.e., it scales as 1/t. In the 
long-time limit, t 3> 1/A/i, the non-equilibrium zero- 
frequency noise dominates, and the electron correlation 



function decays exponentially with time. This is exactly 
the limit where a non-Gaussian markovian noise should 
also be taken into account. 



A. Gaussian noise 

In the context of the noise detection physic o 24 ' 30 the di- 
mensionless counting variable A in the expression (|2ip for 
the FCS generator plays the role of a coupling constant. 
Typically, it is small, A <C 1, so that the contributions of 
high-order cumulants of current noise to the detector sig- 
nal are negligible^ In contrast, in the physical situation 
that we consider in the present paper A = ±7r, imply- 
ing that the shape of the distribution function may be 
strongly affected by high-order current cumulants. Nev- 
ertheless, it is instructive to first consider Gaussian fluc- 
tuations, simply truncating the cumulant expansion at 
second order in A. In this case the correlation function 
(|2T)]) may be evaluated exactly. The are many reasons 
for starting the analysis from considering an example of 
a Gaussian noise: First of all, in equilibrium the current 
fluctuations in a chiral ID system are always Gaussian. 
Second, as we show in the Appendix [Cl the dispersion 
of the charged and dipole modes leads to a suppression 
of higher order cumulants at large distances L. Finally, 
on the Gaussian level it is more easy to investigate and 
compare contributions of zero-point fluctuations and of 
non-equilibrium noise to the electron correlation func- 
tion. 

Thus, expanding the logarithm of the right hand side 
of the Eq. (|13p to second order in A and accounting for 
the Eq. (Ti"4]) . we obtain 



log[xa(A,t)] =i\(j a )t- X 2 J a {t). 



(22) 



Here the Gaussian contribution of current fluctuations 
Sj a (t) = j a (t) — (j a ) is given by the following integral 



J a (t) 



1 

2^ 



dioS a {uj) 

LU 2 + T\ 2 



(l-e-***), ?7^0, (23) 



where the non-symmetrized noise power spectrum is de- 
fined as 



S a {uj) = / dte iut (Sj a (t)6j a (0)). 



(24) 



In what follows, we apply this result for the evaluation 
of the electron correlation function in the case of equilib- 
rium boundary conditions and in the case of a Gaussian 
noise far away from equilibrium. 



1. Equilibrium boundary conditions 

One may propose the following simple test of the non- 
equilibrium bosonization method: Let us consider an 
infinite QH edge. In equilibrium, the charge densities 
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and edge currents exhibit thermal fluctuations. This is 
the case, in particular, for the currents j a through the 
cross-section x = 0, which are considered to be bound- 
ary conditions for the field 4> a in our theory. Therefore, 
one may evaluate the electron correlation function using 
these boundary conditions and compare it with the result 
of the standard equilibrium bosonization technique^ ap- 
plied to a chiral ID system^ 

In equilibrium, (j a ) — 0. The current noise 
power spectrum is given by the fluctuation-dissipation 
relation 31 



S a (u)= / dte wt {j a (t)j a (0)) = 



1 



2tt 1 - e~P u ' 



Substituting this expression into the equation 
obtains 



l0g[ X a(A,i)] 



"4tt 2 



dui 1 



oj 1 



(25) 



one 



(26) 



This integral may be evaluated by expanding the inte- 
grand in Boltzmann factors e ± ^" and integrating each 
term. Substituting the result (for A = n) into Eq. (|2"T|) . 
we arrive at the following expression for the electron cor- 
relation function in the case of equilibrium boundary con- 
ditions: 



K(t) 



sinh(7rt/^) : 



(27) 



which is, in fact, the equilibrium fermionic correlation 
function. The straightforward calculations of the integral 
(lllaj) gives, naturally, the equilibrium distribution func- 
tion, /i(e) = l/(l + e ,9e ) = /f(c) 5 where we have fixed the 
normalization constant, as explained above. Thus for chi- 
ral, interacting quasi-lD systems with a linear spectrum 
equilibrium bosons also implies equilibrium distribution 
of fcrmions. 

It is instructive to compare this result with the known 
expression for the electron correlation function at v = 2, 
found earlier in Ref. [l8[ with the help of the standard 
bosonization technique: 



y-ut^ 



v(3/tt 



u/3/tt 



-i/2 



(28) 

For x — y, details of the interaction leading to wave- 
packet splitting (see Fig. [2]) vanish, and one obtains the 
expression (|2"T1) . thus validating our approach. Moreover, 
the free-fermionic character of the correlation function at 
coincident points (1271) justifies the assumption underlying 
the non-equilibrium bosonization procedure that the FCS 
generators ([TB")) may be taken as for free electrons. 



2. Gaussian noise away from equilibrium 

For a QPC far away from equilibrium, j3A.fi 1, one 
may simply set the temperature to zero. Straightforward 



calculations based on the scattering theory 3 ^ give the fol- 
lowing result for the spectral density of noise (pM|) of a 
QPC: 



S a (u>) = S q (uj) + R a T a S n (u), 



(29) 



where T a = 1 — R a is the transparency of a QPC 
(i.e., the average occupation in the channel a), S q (uj) — 
(1/2tt)oj9(uj) is the quantum, ground-state spectral func- 
tion, and S n (w) = J2± -Sq( w i A/i) — 2Sq(ui), is the non- 
equilibrium contribution (see Fig. Note, that the 
noise power (|29[) differs from the one for a non-chiral 



case 
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'S n 


A.jx 




2tt / 





-A/i 







FIG. 3: Two contributions to the spectral density of noise 
(|29[) . Left panel: Quantum contribution S q (uj) generated 
by the incoming Fermi sea. This contribution vanishes at 
low frequencies S l q (0) = 0, but dominates the behavior of 
the correlator (|23p at short times, tAfi <C 1. Right panel: 
Non-equilibrium contribution S n (u}) dominates at long times 
tAfi ^> 1, i.e., in markovian limit. 

Evaluating the integral (f2"B"| , we arrive at the following 
expression 



J a (t) = (l/4^ 2 )[logt + 2R a T a f(Afit)} 



(30) 



where the logarithm of time originates from the quantum 
contribution Sq, and the dimensionless function /(A/it), 
describing non-equilibrium noise, is given by the integral 



f(Afit) = [ds 
Ja 



[1 — cos(A^s)]. 



(31) 



This function has a quadratic behavior f(Afit) ~ 
(A fit) 2 /4 at short times A fit <C 1, while in the long-time 
(markovian) limit, A fit 3> 1, the dominant contribution 
to this function is linear in time: f(Afit) ~ (it /2)\Afit\. 

For the purpose of further analysis we need the electron 
correlation function in the long-time limit. Taking into 
account that (j a ) = AfiT a /2Tr, we find the cumulant 
generating function 



log[Xo 



2n 
A 

2^ 



AfiT a t 



(log* - nR a T a \Afit\), AfjLt > 1(32) 



Finally, substituting this result into the equation (|21[) . 
and setting T\ = T and T2 = 1 according to the situa- 
tion shown in Fig. [TJ we obtain the electron correlation 
function in the long-time limit: 



K{t) oc t 



-1 iA/iTt- 



-rfTA M | t |/2 ; A/Xt>>1 . 



(33) 
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Note, that the expression (|33| contains the quantum con- 
tribution in the form of a single pole, as for free fermions, 
as well as the non-equilibrium contribution in the form of 
an exponential envelop, those width is determined by the 
noise power at zero frequency, 5i(0) = RTA/j,/2tt. The 
phase shift of the correlator is determined by the "aver- 
age" voltage bias (A/t) = A/iT of the incoming stream 
of electrons, diluted by the QPC. In the next section 
we show that this mean- field like effect of the dilution is 
strongly modified by a non-Gaussian component of noise. 

B. Non-Gaussian markovian noise 

Here we consider non-Gaussian noise and show that 
the contribution of high-order cumulants of current to the 
correlation function is not small. Note, that the quantum 
ground state part of the current noise, S q , that dominates 
at short times, is pure Gaussian. Therefore, the denom- 
inator in the expression (1331) remains unchanged. In the 
long time, markovian limit, the dominant contribution to 
the FCS generator comes from the non-equilibrium part 
of noise, which, e.g., is described by S n in Gaussian case. 
For a QPC, the markovian FCS generator is given by the 
well known expression 24 for a Binomial process: 

Xi(\,t) = (R + Te* x ) N , (34) 

where N = Ant/2ir is the total number of electrons that 
contribute to noise. Applying the analytical continuation 
A — > 7r, we obtain 

log[xi(7r,t)] =^-[]og\T-R\+m6(T-Rj\. (35) 

Substituting this expression to the correlation function 
(f2~Tj) . we arrive at the result 

K(t) OC r V<>(T-R)A^+log iT-illA^IA (gg) 

where the imaginary part of the exponent determines the 
effective voltage bias, while the real part is responsible 
for dephasing. 

Interestingly, at the point T = 1/2 the dephasing rate 
is divergent, and the effective voltage bias drops to zero 
abruptly for T < 1/2. It has been predicted in Ref. 
Q that this behavior may lead to a phase transition in 
the visibility of Aharonov-Bohm oscillations in electronic 
Mach-Zehnder interferometers. We will argue below that 
no sharp transition arises in the electron distribution 
function. However, it leads to its strong asymmetry with 
respect to the average voltage bias (A/i) = TA/j, of the 
outer channel. 



IV. ELECTRON DISTRIBUTION FUNCTION 



analyze the electronic distribution function. We start by 
noting that the experiments [1.2] are done in a particular 
regime of strong partitioning T w 0.5 at the QPC inject- 
ing current to the channel a = 1. This detail, which seem 
to be irrelevant from the first glance, is in fact of crucial 
importance. Indeed, as it follows from the expressions 
(f33|) and (f36|) , the main contribution to the integral (jllaj) 
for the correlation function comes from times t of the or- 
der of the correlation time r c ~ 1/Afi, where our results 
based on the markovian noise approximation are, strictly 
speaking, not valid. However, the numerical calculations 
show that the non-equilibrium distribution in this regime 
is very close to the equilibrium one. Therefore, the ac- 
tual equilibration of electrons, which is reported in the 
experiment [l[ to occur at distances v/ Afi, may in fact 
take place at even longer distances L cq 3> v/ A/i due to 
an unknown mechanism (not considered here). 

Indeed, if the chiral LL model considered in our pa- 
per is valid, then neither the strong interaction between 
electrons of two edges taken alone, nor the weak disper- 
sion of plasmons resulting from a long-range character of 
Coulomb interaction may lead to the equilibration, be- 
cause the systems remain integrable. Thus it seems to be 
reasonable to assume that the equilibration length L c(i 
may indeed be quite long. Therefore, in order to explore 
the physics of collective charge excitations at intermedi- 
ate distances we propose to consider a regime of weak 
injection at the QPC: T -C 1. Firstly, we note that in 
this case our results (|3"3"|) and (|3T>)) may indeed be used 
to evaluate the electron distribution function, because 
the main contribution to the integral (|lla|) arises from 
markovian time scales. Secondly, and more importantly, 
in this regime the electron distribution function acquires 
a strongly non-equilibrium form and the width of the or- 
der of TA/x, which plays a role of the new energy scale. 
Moreover, the advantage of the weak injection regime is 
that it allows to investigate a non-trivial evolution of the 
distribution function, which arises due to bosonic, col- 
lective character of excitations and goes via several well 
distinguishable steps. 

A. Short distances 

At distances of the order of the energy exchange length 
scale 

L cx = v/Afi (37) 

the initial double step distribution function is strongly 
perturbed by the interaction between channels. As we 
argued in Sec. Ill Bl at distances L ^> L ox the charged and 
dipole modes split and make independent contributions 
to the electron correlation function. Therefore, we may 
rely on the result (|36[) . Applying the limit T <C 1 to this 
expression and evaluating the Fourier transform, we find: 



In this section we use the results (|27]) . ([33)) and ([36]) 
for the correlation function of electrons to evaluate and 
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Here, the missing prefactor in the correlation function has 
been fixed by the requirement that /(e) = 1 at e — > — oo. 
Thus, we conclude that energy derivative of the distribu- 
tion function acquires a narrow Lorentzian peak, which 
is shifted with respect to the average bias (A/i) = TA/i 
and centered at e = 0. The last effect is a unique sig- 
nature of the non-Gaussian character of noise. Because 
of the electron-hole symmetry of the Binomial process, 
in the limit J?« 1 the Lorentzian peak obviously has a 
width T ng = 2i?A/j/7r and centered at e = A/i. 

We stress again, that the result (|3"5)) holds only for 
small enough energies close to the Fermi level, namely, 
for |e| < A/i, where the main contribution arises from the 
noise in markovian limit. In fact, the result (|38p fails at 
large energies in somewhat non-trivial way. Namely, it is 
easy to see that any electron distribution function has to 
satisfy the sum rule 



(A/i) = e + / def(e) 



deedf(e)/de, (39) 



where eo is the cutoff well below the Fermi level, and 
the "average" bias (A/i) = TA/i in the case of linear 
dispersion of plasmons. This sum rule simply expresses 
the requirement of the conservation of the charge current 
and implies certain amount of asymmetry in the distri- 
bution function. In the present case, such an asymmetry 
arises in the power-law tails of the function —df(e)/de 
and originates from quantum non-equilibrium noise. It 
can be seen in Fig. [H where the results of numerical cal- 
culations are shown. 

Moreover, at energies of the order of A/i the power- 
law behavior of the function (f38|) has to have a cut-off, 
because the QPC does not provide energy much larger 
than the voltage bias. Quantitatively, this follows from 
the conservation of the energy. We demonstrate below 
that for the system with linear dispersion of plasmons, 
the heat flux in the outer channel can be written entirely 
in terms of the single-electron distribution function (in 
unites e = h = 1), 



7 m = (1/2tt) J de< 



[/(e)-0((A/i)-e)] 



(40) 



as in the case of free electrons. We use the subscript 
"m" in order to emphasize the fact that it is this quan- 
tity that has been measured in the experiment [2|. In 
Sec. [V] we show that at distances L ^> i cx the total heat 
flux injected at a QPC splits equally between two edge 
channels, therefore integrating Eq. (|4T)|) by parts and sub- 
stituting the heat flux for a double-step distribution, we 
obtain 



/n 



(TA/i) 2 

47T 



1 

47T 



dee 2 



df(e) _ TR(Atf 



dc 



Stt 



(41) 



for L ^> T cx . One can see from Eq. (|38|) that indeed the 
power-law behavior has to have a cut-off at |e| ~ A/i. We 
stress, however, that this summation rule is less universal 



than the one given by Eq. (|39|) . because it accounts only 
a single-particle energy of electrons and fails in the case 
of a non-linear spectrum of plasmons, considered in Sec. 
E| in detail. 



B. Intermediate distances 

So far we have considered the case of a linear spectrum 
of plasmons. This is a reasonable assumption, taking into 
account the fact that non-linear corrections in the spec- 
trum of plasmons lead to a nonlinear corrections in the 
Ohmic conductance of a QPC. However, the experiments 
seem to be done in the Ohmic regime. Nevertheless, 
even in the case of a weak non-linearity in the spectrum 
of the the both modes of the sort^ 

kj(uj) — uj/vj + 7jO; 2 sign(a;), v% — u, v-i — v, (42) 

barely seen in the conductance of a QPC, an intermediate 
length scale L g may arise at which high-order cumulants 
of current are suppressed, and the noise becomes effec- 
tively Gaussian. This situation occurs when the wave 
packets of the original width w/(TA/i) overlap. A sim- 
ple estimate using the nonlinear correction (|42|) gives the 
length scale 



T g = l/ 7 (TA/i) 2 , 7 = min( 7 ,). 



(43) 



We support this conclusion by rigorous calculations in 
Appendix [Cj 

The non-linearity in the spectrum is weak, if jvTAfi <C 
1. This implies that L g ^> T cx , and leads to the 
possibility to observe non-Gaussian effects at distances 
L cx <C L <C L s , discussed in the previous section. Ob- 
viously, the same requirement also guarantees that dis- 
persion corrections to the Ohmic conductance of a QPC 
are small. This allows us to neglect corrections to the 
quantum part of the electron correlation function and to 
use the result (|33|) for a Gaussian noise. Substituting this 
result to the equation (jllaj) . we obtain 



am 

de 



( e -(A/i» 2 



r 2 ' 

g 



r g = ^TA/i/2, (44) 



in the case L g <C L <C L cq . One can see, that the width 
of the function (|4"4")l is almost twice as large compared to 
the one the function (|38p. Moreover, the function (|4"4")) 
satisfies the sum rule (|39]) . Therefore, we do not expect 
any asymmetry in the high-energy tails of this function, 
in contrast to the situation with the non-Gaussian noise. 
The comparison of distribution functions in these two 
regimes is shown in Fig. [4] 

So far we have considered a situation where both 
charged and dipole mode are dispersive. If for some rea- 
son the dispersion of one of the modes is negligible, then 
higher order cumulants are suppressed only by a factor of 
two. The derivative of the electron distribution function 
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FIG. 4: Energy derivative of the electron distribution func- 
tion, —df/de, is shown for different distances L from the 
QPC injecting current. The transparency of the QPC is set 
to T = 0.05 and voltage bias is A/j, — 40 fiV. Black line: 
—df(e)/de for short distances L cx < L < L s , so that the 
noise is Gaussian (|38[) . Red line: —df(e)/de for intermediate 
distances L g < L < L oq , where the noise is Gaussian (|44|l . 
Blue line: The derivative of the Fermi distribution function 
at the temperature given by Eq. (|47fl . The dashed line is a 
guide for eyes at the energy equal to the effective voltage bias 
(A/x) = TAfi = 2fiV. Inset: The same distribution functions 
are shown in the integrated form. They are shifted vertically 
by 0.2 intervals for clarity. 

in this situation is given by the Lorentzian 

dm = (r ng + r g )/27r 

de (e-<A M )/2)2 + (rV + r g )74 1 ' 

centered at (A/z)/2 = A/iT/2 with the width (T ng + 
r g )/2 = (1/tt + 7r/4)A/iT. This is because one mode 
brings only Gaussian component of the markovian noise, 
while the other one brings full non-Gaussian noise. 



C. Long distances 

Next, we consider the distribution function at long dis- 
tances, L L cq , after the equilibration takes place. The 
temperature of the eventual equilibrium distribution may 
be found from the conservation of energy. In the next sec- 
tion we show that the heat flux produced at QPC splits 
equally between two edge states. In the situation of linear 
dispersion the distribution function acquires the form 

= 1 + e(e -lAM))/r oq - ( 46 ) 

The possibility of such equilibration process is suggested 
by the fact that the equilibrium distribution of bosons 
implies the equilibrium distribution for electrons, as has 
been shown in Sec. IIII A 11 Obviously, the distribution 
(l4T>)) satisfies the sum rule (|59")l . while the energy conser- 
vation condition (1411) may now be used in order to find 



the effective temperature: 

T cq = V3T/2^ 2 A M , (47) 

where we have used T<1. 

We conclude that the width of the equilibrium distri- 
bution scales as vT, in contrast to the case of a non- 
equilibrium distribution at shorter distances from the 
current source, where it scales linear in T. Therefore, 
if T is small, an equilibrium and non-equilibrium dis- 
tributions may easily be distinguished, as illustrated in 
Fig. In the situation where the dispersion can not be 
neglected, the equilibrium distribution of fermions is not 
given by the Fermi function (|46l) . This situation deserves 
a separate consideration, which is provided in the next 
section. 



V. MEASURED AND TOTAL HEAT FLUXES 

We have seen that in the case of weakly dispersive 
plasmons, jvTA/j, <C 1, the non-linearity in the spec- 
trum leads to the suppression of high-order cumulants of 
current noise at relatively long distances, which strongly 
affects the distribution function. On the other hand, 
the direct contribution of the non-linear correction in the 
spectrum to local physical quantities, such as the QPC 
conductance and the heat flux, is small and has been so 
far neglected. Nevertheless, it may manifest itself exper- 
imentally in a quite remarkable way. In this section we 
show that the non-linearity in the the plasmon spectrum 
contributes differently to the measured heat flux ([4TJ]) and 
to the actual heat flux expected from the simple evalu- 
ation of the Joule heat. As we demonstrate below, this 
may, under certain circumstances, explain the missing 
energy paradox in the experiment [2j . 

We start by noting that the measured flux (|4TJ|) at the 
distance L form the QPC may be expressed entirely in 
terms of the excess noise spectrum S a (cj) = S a (uj) — 
(1/2tt)ui9(ui) of edge channels right after the QPC, where 
S a (u>) is defined in (|2~4"]l . Namely, in Appendix [B] wc 
derive the following result: 

1 f°° 

I m (L) = - / dw{Si(w)[l + cos(AfcL)] 

J — OO 

+ S 2 (w)[l -cos(AfcL)]}, (48) 

where Afe = ki(u>) — k-z{uS), and kj{—w) = —kj(u>). Im- 
portantly, this result holds for an arbitrary non-linear 
spectrum kj(uj) of the charged and dipole modes, and for 
a non-Gaussian noise in general, i.e., high-order cumu- 
lants simply do not contribute. 

One can easily find two important limits of Eq. (|48l) : 
for L = we immediately obtain an expected result 

lm(0) = ~ / dwSi(w), (49) 
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while at L — > oo the cosine in (l48l) acquires fast oscilla- 
tions, and we get 

Im(oo) = l [ dw[Si(w)+Sa(w)]. (50) 

To be more precise, this happens at L > L cx = w/A/z. 
At zero temperature §2 vanishes, and the single-electron 
heat flux 7 m , created at the QPC, splits equally between 
edge channels: / m (oo) = I m (0)/2. Note also, that in 
the case of linear dispersion § Q = T a R a S n , where S n is 
shown in Fig. [3l 

In the next step, we rewrite the same measured flux in 
terms of the plasmon distributions rij(k) — (a^(k)dj(k)), 
see Appendix [Bl 

1 f dk 

J -(°°) = — £ / T uj](k)n 3 (k) + I q , (51) 

j 

= ^E/fHW-M)^ (52) 
i J 

where Vj — dujj/dk are the plasmon speeds at k = 0. 
The term I q is the contribution to the measured flux 
from quantum smearing of the zero-temperature electron 
distribution function / (e) close to the Fermi level, which 
originates from a non-linear dispersion of plasmons. 

Here an important remark is in order. The integral 
(|52f may diverge at large k and has to be cut off at the 
upper limit. This is because there is no guarantee of the 
free-fermionic behavior of the correlator K(t) at short 
times and of the zero-temperature electron distribution 
function /o(e) at large energies. Thus, the integral (|4"0")) 
has to be also cut off, which is what in fact is done in 
experiment. In contrast, the spectrum of plasmons is 
linear at small k, and thus the distribution function /o(e) 
has a free-fermionic behavior close to the Fermi level. 
Our definition of I q corresponds to the normalization of 
/o(e) to have a discontinuity of the value — 1 at e = (A/i). 
The experimentally measured I q may differ from the one 
defined in (|52[) by a constant, which is, on the other hand, 
independent on the voltage bias A/x. 

Next, we note that the actual total heat flux in the 
case of a non-linear dispersion of plasmons acquires the 
completely different formal 

/h = ^E jdk^^n^k), (53) 
3 

and thus in general I m =/= /h/2, contrary to what has 
been assumed in the experiment Q. This may explain 
the missing energy paradox. Indeed, assuming the low w 
spectrum of the general form 

k 3 = u/vj +i 3 u l i, j = l,2, (54) 

where jj are small, and equilibration of plasmons at L — > 
00, i.e., nj(k) = n B (uj/T eq ) = l/[exp(w J -/r eq ) - 1], we 
obtain the missing heat flux as 

Im - A - 4/2 = J2 c n^ 3 T%+\ (55) 

j=l,2 



where the constants Cj = (1/An) J dzz lj+1 riB{z) are of 
the order of 1. This result implies that experimentally, 
the missing heat flux may be found investigating its bias 
dependence and the spectrum of plasmons. 




FIG. 5: Typical spectrum of charged plasmon in the case of 
the Coulomb interaction screened at distances d 1/fc. This 
spectrum is concave, i.e., duj/dk < ui{k)/k. 

Let us consider an example where the dispersion of 
charged plasmon at small k arises from the screened 
Coulomb interaction^ 

Lj/k = 2[JC (ka)-IC (kd)] 

= 2\n(d/a) - (l/2)(fcd) 2 tn(2/fcd), (56) 

where a is the high-energy cutoff, d is the distance to 
the gate, such as ka <§C kd -C 1, and /Co is the MacDon- 
ald function. The low- A: asymptotics of this spectrum is 
illustrated in Fig. [5j One can see that the spectrum is 
concave, so that in this case the measured heat flux (|5"Tj) 
is larger than the half of the actual heat flux ([53)1 . In 
addition, the effect is weak, because kd ~ 0.1 in the ex- 
periment 0. Thus, the dispersion of the Coulomb inter- 
action potential alone is not able to explain the missing 
flux paradox. Various mechanisms of convex dispersion 
are still possible and will be investigated elsewhere. 

VI. CONCLUSIONS 

Earlier theoretical works on quantum Hall edge states 
at integer filling factors may be divided into two groups: 
fermion based and boson based theories. Recent in- 
terference experiments suggest that the boson approach 
might be more appropriate for the description of the edge 
physics. However, both groups of theories give the same 
predictions for the local physical quantities at equilib- 
rium. Moreover, the first theoretical works based on 
fermion^ and boson 4 approaches and addressing the non- 
equilibrium local measurements, have not been able to 
make qualitatively distinct predictions. In this paper 
we show that it is nevertheless possible to test and dif- 
ferentiate between two approaches with the local non- 
equilibrium measurements. 

We address recent experiments [l||2( with quantum Hall 
edge states at filling factor 2, where an energy relax- 
ation process has been investigated by creating a non- 
equilibrium state at the edge with the help of a QPC and 
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reading out the electron distribution downstream using a 
quantum dot. We use the non-equilibrium bosonization 
approach 5 in order to describe the gradual relaxation of 
initially non-equilibrium, double-step electron distribu- 
tion to its equilibrium form. In the framework of this ap- 
proach the non-equilibrium initial state is encoded in the 
boundary conditions for the equations of motions that 
depend on interactions. We show that the electrons ex- 
cite two plasmons: fast charged and slow dipole mode. 
Thus the electron correlation function (f2~Tj) is expressed 
in terms of the four contributions, each having form of 
the generator of FCS of free electrons with the coupling 
constant A = n. Evaluating the Fourier transform of this 
function, we find the electron distribution function. 

Before reaching eventual equilibrium form, the distri- 
bution function evolves via several steps, where its energy 
derivative acquires a Lorentzian shape: 



dm 

de 



T/ir 



< 



Am, 



(57) 



Here the width T and centering eo take different values 
in different regimes. Each of the regimes, summarized 
below and illustrated in Fig. [5J has its own dominant 
process: 

(i) First, after tunneling through the QPC, electrons 
excite plasmons, which then split in two eigenmodes: one 
is charged fast mode with the speed u, and the other is 
slow dipole mode with the speed v. This process takes 
place at distances L cx = v/Afj,, where Afi is the voltage 
bias across the QPC. In this regime Eq. (|21l) applies, 
which eventually leads to the the distribution ([57} with 
the width T = T ng — 2A/iT/7r, centered at eo = 0. 

(ii) Next, a weak dispersion of plasmons, e.g., of the 
form k = eo / v + juj 2 sign(uj) , leads to broadening of wave- 
packets of the energy width e and to their overlap. This 
process takes place at distances L ^> 1/je 2 . As a result, 
high-order cumulants of the current injected at the QPC 
are suppressed at distances L 3> L g — l/j(TAfi) 2 , the 
noise becomes Gaussian, and the derivative of the elec- 
tron distribution function acquires the shape (|57[) with 
the width r = T g = irAfiT/2, centered at e = AfiT. 

(iii) A situation is possible, where the dispersion of one 
mode, most likely of the charged mode, is much stronger 
than the dispersion of the second mode, i.e., 71 ^72- In 
this case, the previously described regime splits in two 
separate regimes. First, at distances L — l/ji(TA[i) 2 
the contribution of the charged mode to high-order cumu- 
lants of noise become suppressed, which leads to the dis- 
tribution (|57l) with the parameters T = (r ng + T g )/2 and 
e = A^iT/2. Then, at longer distances L = 1/7 2 (TA^) 2 
the noise becomes fully Gaussian. 

(iv) The interaction may lead to broadening of the 
wave-packets, but they do not decay, which implies that 
the interaction alone does not lead to the equilibration. 
This means that a different, weaker process may lead to 
the equilibration at distances L eq much longer than the 
discussed above length scales. In the tunneling regime, 
T <C 1 the width of the eventual equilibrium distribution 



scales as y/T, in contrast to the above regimes, where it 
scales as T. Thus, to observe the described variety of 
regimes, we propose to perform measurements at large 
voltage biases and low transparencies of the QPC uti- 
lized to inject electrons. 

Finally, we suggest a possible explanation of the para- 
dox of missing heat flux in the experiment [2j . So far we 
have summarized the effects of weak dispersion, which 
lead to appearance of intermediate length scales. We 
have found that in the case of a strongly nonlinear dis- 
persion of plasmons, the measured heat flux I m in the 
outmost edge channel, experimentally determined with 
the procedure described by Eq. (j40j) . is different from the 
actual heat flux per channel Ih/2, defined by Eq. ([53]). 
The screened Coulomb interaction leads to a rather weak 
dispersion of the charged plasmon, and the effect is of the 
opposite sign, because the spectrum in this case is con- 
cave. Nevertheless, other mechanisms of the convex dis- 
persion are possible. They will be considered elsewhere. 



-*~ \. ■* oc v'/'A/i 




■j/An 1/7A/X 2 1/7TV L, 



eq 



FIG. 6: Different length scales for energy relaxation processes 
and corresponding distribution functions in each regime are 
schematically shown. Red curve: The initial double-step dis- 
tribution function. Black curve: At distances L 3> Afi/v the 
distribution function is strongly asymmetric with respect to 
the "average" bias (A/it) = TAfi. Green curve: At distances 
L 3> l/7(TA/i) 2 the distribution function is a Lorentzian 
with the width that scales as TA/i. Blue curve: The final 
equilibrium Fermi function at large distances. For small trans- 
parencies its width scales as VtA/i. 
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Appendix A: Solution of equations of motion 



In this appendix we solve the equations of motion (|12a|) 
with the boundary conditions (| 12h>[) in the case of a po- 
tential U a p(x — y) of a finite range, where the plasmon 
spectrum is non-linear. For doing this, we first write 
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down the normal mode expansion for the edge boson 
fields: 



+ 



ikx I t 



-ikx] 



(Al) 



We consider zero modes to be classical variable, because 
the commutator [ir a , (p a ] — i/W vanishes in the thermo- 
dynamic limit W — > 0. Then, we rewrite the operators 
a ak in the new basis dj k , which diagonalizes the edge 
Hamiltonian ([8]): 

a Jk (t) = Oj^W, (A2) 

where j = 1,2, and u>j(k) is the dispersion of the jth 
mode. In the case, where the in-channel interaction 
strength is approximately equal to the intra-channel one, 
U a p(x — y) fa U(x ~ y), and the interaction is strong, 
U(k) = J dxe lkx U(x) ^> 27TUf, the transformation to the 
new basis is simple and universal: 



aifc(t) 



= (a lk e 



+ ~a 2k e~ l ^ (k)t ), (A3a) 



a 2k (t) = -J= {aike-^W - 5 M e^ 2 ^). (A3b) 

Then, wi(fc) = k[vp + U(k)/ir] for the charged plasmon, 
and W2{k) — vk for the dipole mode, where v pa vf- 

In the next step, we use the boundary conditions (|12b[) 
to connect the current operators to the operators (|A2|) : 



i 9wi irk , . , , . , /ha \ 

aik = ~ 1 ^kyw biiui)+j2{uJi)i (A4a) 

Z / 7r/c 



" 2fe = ^ ~dk\l w Mu}2) - j2 ^> (A4b) 

where j a (uj) — J dte ttJJt j a (t) , and we have used the ob- 
vious relation j a (— lo) = j^u). Finally, substituting re- 
lations (|A4[) into Eq. (|A3ap and then to the expansion 
(lAip . we find the solution of the equations of motion for 
the boson fields. In particular, 



fa(x,t) = -2n(j 1 )t + 



+ n(u)(e tklX - e^ x )}e 



(A5) 



where we set kj(-uj) = —kj{uj). In addition, we have 
omitted the contribution of the zero mode 7Ti, because 
we need local correlators, and replaced the zero mode 
ipi(t) by its expectation value. 



Appendix B: Evaluation of the measured heat flux 

The experimentally found heat flux, defined in Eq. 
(1401) . may be written in time representation as 



I m = -id* 



K(t)- 



,i{An)t 
2-nit 



(Bl) 



Then, we may use Eq. (jllbp and results of the Appendix 
lAlin order to evaluate the correlation function K(t). The 
difficulty of finding this function is related to the fact 
that according to Eq. (|A5I) the correlators of the oper- 
ator fa(t) are determined by the currents j a {t), which 
are in general non-Gaussian. However, in the present 
case we need to take the limit t — >• in Eq. (|B1I) . The 
high-order cumulants of currents originate from a non- 
equilibrium process and are suppressed at short times 
tA[i <C 1. Therefore, we are allowed to evaluate K(t) in 
Gaussian approximation. 

This may be done by expanding the right hand side of 
the Eq. (jllbp to second order in fa: 

lnif(i) = ([fa(t) - fa(0)}M0)) + 27TtO'i>t, (B2) 

where the averaging is over the fluctuations of the cur- 
rents j a . We then use the Eq. (|A5|) and the stationary 
correlators of the currents 

(ja(ui)jp(u2)) = 2TrS a (uj 1 )S al3 S(uj 1 + u 2 ) (B3) 

to present the electron correlation function in the follow- 
ing form 

hxK = lnK n + 2ni(j 1 )t - Int. (B4) 
Here the fluctuation contribution reads 

oo 

\nK n = TT / ^(e-' i " t -l){§i(w)[l + cos(Afc J L)] 



+ § 2 (w)[l-cos(AfcL)]}, (B5) 

where we have introduced the excess noise spectral densi- 
ties S a (w) = S a (u)) — v9(v)/2w and Ak = fci(w) — k 2 {oj). 

It is easy to see that for a non- vanishing contribution 
to I m , given by the expression (1B1|) . we need to expand 
lnifn to second order in t. Note, however, that linear in 
t term in this expansion adds to the corresponding term 
in Eq. (|B4I) to give 2m(Afi)t. This follows directly from 
the definition ([39)) of the "average" bias, which in time 
representation may be written as (A/i) = 2irdt[tK(t)]t=o- 
Therefore, only t 2 term in hiK n contributes to I m , and 
we obtain 



1 

47T 



d 2 \nK n (t)\ t=0 . 



(B6) 



Finally, we use Eq. (|B5[) and obtain the result (|48[) . 

Next, we wish to rewrite the measured flux I m in terms 
of the plasmon distributions rij(k) — (at(fc)aj(fc)), j — 
1, 2. For doing so, we now use Eqs. (jAll) and (IA3|) . repeat 
the steps that lead to (|B5|) . and take the limit of L ^> L cx . 
The result may be presented in the following form: 

oo 

/dk 
— rij(k)[l — cos(ujjt)] 

3 o 

oo 

+ IY.Jt (e-^'-e"™^*), (B7) 



J o 
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where the last term is the quantum contribution due to 
the non-linear plasmon spectrum, and Vj — duij/dk are 
the plasmon speeds at k = 0. Substituting expressions 
(|B7| into Eq. (|B6)l . we obtain the final result ([5T]) and 
(l52l) for the measured heat flux. 



Appendix C: Suppression of higher order cumulants 



One can easily see that Y^i^ji^i) — 7j U} fsign-(wi), 
because the integrals in (|CJ3|) are limited to J^w; = 0- 
For the second cumulant this implies that kj(uii) + 
kj(u)2) = 0, i.e., the dispersion correction cancels too. 
Therefore the second cumulant is not suppressed at long 
distances. Below we consider high-order cumulants. Us- 
ing the identity 2tt8(uji + . . .+oj n ) = J dr exp[i(u;i + . . . + 
w„)t] we can write 



Here we show that a weak dispersion of plasmon modes 
leads to the suppression of the contribution of higher 
order cumulants of current j a to the electron correla- 
tion function at long distances L from the source of 
currents. We demonstrate this using an example of 
a weakly dispersive spectrum of plasmons in the form 
kj = cu/vj + 7j-u; 2 sign(u;), j = 1,2. Since we are inter- 
ested in the behavior of the electron distribution function 
close to the Fermi level, we need to know a long-time 
asymptotics of the electron correlation function K(t). 
Therefore, the contributing currents j a can be consid- 
ered Markovian processes and fields <fi a can be treated as 
classical variables. 

Let us consider the nth cumulant: 

M^(L,t) = (([ML,t)-MLM n ))- (CI) 

According to Eq. (|A5|) . at large distances L ^ L cx = 
Vj/ Afi and long times tAfi 3> 1 it may be written as 



MW(I,f) = ^lW(L,f), 



(C2) 



oo 

MtHL,t) = S^Yl J driF^LT, (C4) 



where we have introduced the integrals 

p * f I „— iuit _ i \ iujT+i~tjLu 2 sigTi(u]) 

J ~ 2 J oj [C )S 



(C5) 



At large distances L^j ^ t 2 the contribution to the inte- 
grals Fj comes from small w, where one can approximate 
e -iuit _ j ~ —iut. Therefore, Eq. (IC5[) can be further 
simplified: 

Fj = (t/2) / duje luJT+l ^ Luj2sien( - UJ ' ) oc _^ e ±4T2/4 ^ L . 

(C6) 

Substituting this result into Eq. (|C4[) and then to (|C2I) . 
we find that 

MW(L,t)« tJ2St } Y, n>2 ' ( C7 ) 



a j=l,2 



where: 



M<j?\L,t) = 2irS^ j f[ ^i(i/2)(e- Mlt - 1) 
J i=i U>1 

x 6(0)! + ... + uj n ) \e 1 ^ kl{ui)L + e*£i k2(uJl)L (C3) 

and = ({ja}}- Here we have neglected the cross terms 
containing fast oscillating functions. These terms have 
the same origin as fast oscillating terms in (|4"5)l and van- 
ish at distances L 3> L cx . Dropping those terms is also 
equivalent to neglecting in (|2"U)) correlations of charges 
taken at different times t u and t v . Finally, we note that 
in our particular case, where the QPC is connected to 

(n) 

the outmost edge channel only, 5^ = for n > 2. 



where, we recall, the sum is over the plasmon eigenmodc 
number j and over the channel number a. 

We note, that at large distances L the cumulants 
M^(L, t) are suppressed by the dimensionless small pa- 
rameter t 2 /jjL <C 1. In our case s[ n ^ ~ TAfi and 



S. 



(n) _ 



0, therefore the contribution of high-order cumu- 



lants to the correlator K(t) may be neglected at distances 
larger than L g — l/jj(TA[i) 2 , and noise may be consid- 
ered Gaussian. Obviously, if only one plasmon mode of 
two is dispersive, e.g., 72 = 0, then at distances L L g 
the cumulant (|C2[) is suppressed by the factor of 2. One 
can interpret the result (|C7|) as the renormalization of 
the effective coupling constant A in (fl~3)) . which is caused 
by spreading of a plasmon wave packets due to the dis- 
persion. 
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